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Abstract 

CUORE is a 1 ton scale cryogenic experiment aiming at the measurement of 
the Majorana mass of the electron neutrino. The detector is an array of 988 
Te02 bolometers used for a calorimetric detection of the two electrons emitted 
in the /3/3(Oz^) of -'^'^^Te. The sensitivity of the experiment to the lowest Majorana 
mass is determined by the rate of background events that can mimic a /3/3(Oz/) . In 
this paper we investigate the contribution of external sources i.e. environmental 
gammas, neutrons and cosmic ray muons to the CUORE background and show 
that the shielding setup designed for CUORE guarantees a reduction of this 
external background down to a level < 10~^ counts/kg keVyr at the Q-value, as 
required by the physical goal of the experiment. 

Keywords: Neutrinoless Double Beta Decay, Cosmic Ray Background, 
CUORE 



1. Introduction 

The sensitivity of an experiment looking for very rare events such as /3/3(Oj^) can 
be spoiled by background counts in the region of interest (ROI). The most dan- 
gerous background sources for such experiments are due to cosmogenic activa- 
tion and intrinsic contamination of the detector and setup materials, and to 
environmental sources in the experimental hall. In a deep underground labora- 
tory the environmental background can be ascribed to three different sources: 
gammas, neutrons and cosmic ray muons. Gammas and neutrons are due both 
to the natural radioactivity of the rock constituting the laboratory walls and 
to muon interactions in the rock, in the materials surrounding the detector and 
in the detector itself. The reduction of the gamma and neutron background in 
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an underground experiment is generally achieved by means of passive shields 
while a muon identification system is used to tag and reject muon induced 
events. Shields and vetoes are designed to match the specific experimental fea- 
tures, mainly by considering the kind of physical results the experiment aims 
to achieve and considering the various radioactive sources that could spoil the 
experimental sensitivity. 

In this paper we discuss the impact of the above mentioned external sources 
on the sensitivity of the CUORE experiment. In section [5] we describe the ex- 
periment and its physical goals, in section [3] we list the available information 
about the environmental background in the experimental site, while the fol- 
lowing sections are devoted to the description of the Monte Carlo simulation 
code used for the background estimates and of the results obtained in terms of 
contributions to the CUORE background from the different external sources. 



2. The CUORE Project 

The nature of neutrino mass is one of the frontier problems of fundamental 
physics. Neutrinoless Double Beta Decay {Pf3{0i/)) is a powerful tool to in- 
vestigate the nature of the neutrino and in the Major ana hypothesis the mass 
hierarchy [l| and possible extensions of the Standard Model. The /3/3(Oi/) is 
a rare spontaneous nuclear transition where a nucleus (A,Z) decays into an 
(A,Z-|-2) nucleus with the emission of two electrons and no neutrinos, result- 
ing in a monochromatic peak in the sum energy spectrum of the electrons. The 
PPiQv) decay is possible only if the neutrino is a Majorana massive particle. The 
transition width is proportional to the square of |(m,y)|, the effective Majorana 
mass of the electron neutrino. From the f3p{Qv) half-life it is therefore possible 
to infer the correct mass hierarchy and the absolute mass scale of neutrino: 
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where G^'^ is a phase space factor including the couplings, and M^,^ are 
the Fermi and Gamow- Teller nuclear matrix elements, respectively, and 5 a and 
gv are the relative axial vector and vector weak coupling constants, respectively. 

The CUORE (Cryogenic Underground Observatory for Rare Events) exper- 
iment aims at searching for Pj3{Qv) of ^'^"Te with a background < 10"^ 
counts/kg keV yr and an effective Majorana mass sensitivity of a few tens of 
meV, reaching the so called inverse hierarchy region of the neutrino mass spec- 
trum [ij. It is presently under construction in the Hall A of the Gran Sasso 
Underground Laboratory (LNGS) at an average depth of '^3600 m.w.e. Q. 



2.1. The Detector 

The CUORE detector is a system of 988 cryogenic bolometers. The array 
is composed by 19 vertical towers cooled by means of a dilution (■^He/'*He) 
refrigerator to a temperature of ^10 mK. A tower consists of 13 layers of 4 
bolometers each. The bolometer consists of a 5x5x5 cm"^ Te02 cubic crystal. 
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an NTD Ge thermistor, and a heater (a Si resistor doped with As) glued with 
Araldit Rapid Epoxy onto the crystal surface 0]. The dielectric Te02 crystal 
has a very low heat capacity and a very large (27% in mass) natural isotopic 
abundance of the /3^(0j/) candidate ^^°Te . Each cube has a weight of 750 g, 
so the detector is a granular calorimeter of total mass 741 kg, corresponding to 
~200kg of i30Te . 

The cryogenic detecting technique was proposed for the first time by Si- 
mon Q to study nuclear phenomena and then suggested to search for rare 
events by Fiorini and Niinikoski 5]. 

The typical performance of such detectors has been measured in previous ex- 
periments: MIDBD the first experiment operating a large array of bolometers 
to search for /?/3(0^) , CUORICINO 0, a single CUORE tower running for five 
years at LNGS and other tests performed with CUORE prototype detectors. 
They have shown an average energy FWHM resolution of 1-2 keV at 10 keV 
and 5 keV at 3 MeV and an extremely slow time development of the signal 
pulse. The 90-10% average rise time was evaluated to be of about 50 ms, and 
the 90-30% decay time of about 100 ms, with a low tail in its distribution that 
extends even above one second. Typical energy thresholds for double beta decay 
searches are kept between 50 and 300 keV. 

The CUORE cryogenic apparatus is shown in Fig. [TJ The cryostat is made 
of six nested vessels and its base temperature without heat loads is expected 
to be as low as 6mK. Three lead shields are used to protect the detector from 
environmental radioactivity and from contaminations in the building materials. 
A 25 cm thick lead layer outside the OVC shields the detector from radiations 
coming from the bottom and from the sides. An equivalent shielding against 
radiation coming from the top has to be placed inside the cryostat, just above 
the detector. This is a 30 cm thick lead disk with a diameter of about 90 cm. 
Just below it copper disks totaling an additional 8 cm shields are placed. A lead 
ring-shaped shield placed on the Steel flange closes the gap between the lead 
disk above the detector and the outer room temperature shield. An additional 
shielding of detector's sides and bottom is provided by a 6 cm lead layer just 
outside the Steel shield. 

Outside the external lead shield an 18 cm thick polyethylene layer will be 
added in order to thermalize environmental neutrons that will then be absorbed 
by a 2 cm layer of H3BO3 powder contained in the hollow space between the 
lead and the polyethylene itself. On the top and on the bottom of the setup, 
due to construction difficulties, the neutron shield will be made of a single 20 
cm layer made of polyethylene with 5% in boron. 

2.2. Physics goals 

Operated as a single particle detector, the cubic bolometer will be able to 
detect a fii^{Ov) occurring in its bulk by recording the two electrons emitted 
in the decay. The two electrons (and the recoiling nucleus) deposit their whole 
energy in the crystal where the decay takes place with a probability of 86.4% 
(computed with a Monte Carlo simulation). In this case, the bolometer records 
a monochromatic line at the /3f3{Qi^) transition energy (2527.517 ± 0.013 keV 
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Figure 1: 3D view of the cryogenic apparatus. 



for ^-^^Te [1, Q), while the other bolometers of the array are not involved in the 
decay. Background events that can mimic a /3/3(Oi^) decay and hence spoil the 
sensitivity of the experiment are those that correspond to an energy deposition 
similar to that of /3f3{0i'). In CUORE, a large fraction of those events can be 
rejected simply requiring the signal to be a single-hit event as it is expected for 
a /3/3(Oi/) decay. The reduction factor depends on the actual source kind and po- 
sition as shown in Tab.|3l Monte Carlo simulations show that this reduction can 
range from a factor 2 for bulk contaminations of the shields and crystal surface 
contaminations to two orders of magnitude for the cosmic muon contribution. 
The probability of a spurious time coincidence between a background event in 
one crystal of the CUORE array and a PPiQv) event in another crystal is Pcoinc 
= AT-Ndet-/, being Ar=50 ms a reasonable value for the coincidence window, 
corresponding to the signal rise time, Ndet=988 the number of detectors and / 
the average rate. Assuming for / the value measured on CUORICINO detectors 
of 5 mHz we obtain the conservative estimate Pcoinc ^ 25%. This probability 
can be reduced by means of a near channel selection, without loosing in terms of 
background rejection capability. Taking into account only the 18 closest crystals 
the probability of spurious coincidences is reduced to 0.45%. 

The experimental sensitivity is conventionally parametrized in terms of the 
number of /3/3 emitters (N^^), the detection efficiency (e), the live-time (T), the 
energy resolution (F) and the background counting rate {B): 



Soupp = ln2 X N^^ X e X \/ (2) 
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The background counting rate of CUORE is obviously the ultimate param- 
eter that will determine the /3/3(Oz^) sensitivity. CUORICINO has measured a 
background level of 0.18 counts/kg keVyr with a 40.7 kg array of CUORE-like 
detectors. A smaller array consisting of eight 5x5x5 cm^ detectors assembled 
in a two planes CUORE-like tower, after a dedicated surface treatment of the 
copper structures and of the crystals, has been recently operated obtaining a 
background rate of 0.06 counts/kg keVyr . CUORE points toward an even bet- 
ter background level (< 10~^ counts/kg keV yr ) exploiting its granularity, the 
specific design of detector, cryostat and shields and the careful selection and 
preparation of materials. 

3. External background sources 

As already stated, the possibility to study rare events such as /3/3(0:^) is 
strongly influenced by the background in the ROI. There are various sources 
that give rise to spurious counts in this region such as environmental gamma ra- 
dioactivity, cosmic rays, neutrons, radon and contamination of materials which 
detectors and their shielding are made of. The evaluation of the background 
due to cosmogenic activation and intrinsic contamination of the detector com- 
ponents is described in Radon emanation in the experimental hall is not 
an issue for CUORE since all the gamma lines of the^^^i?n have energies below 
the ^^°Te /3/3(Oz^) Q-value. In this paper we focus on the contribution coming 
from external background sources: environmental gamma radioactivity, cosmic 
rays (muons) and neutrons. 

3.1. Gammas 

The gamma spectrum due to natural rock radioactivity at LNGS extends 
up to 3 MeV while above that energy the main contribution to the gamma 
background comes from photons produced by neutron and muon interactions in 
the rocks. Due to their high energies and lower flux (estimated to be 10^ times 
smaller than the flux of the 2.6 MeV gamma-ray of ^os^j rock [llj) these 

photons cannot be measured with a small Ge diode. In Fig. [5] the gamma-ray 
spectrum, as measured by a small Ge diode in the Hall A of LNGS, is shown. 
The reconstructed flux is reported in Tab. [T] [12]. The integral gamma ray flux 
below 3 MeV has been reconstructed to be ~0.73 7/scm^ . 



3.2. Neutrons 

In a deep site the neutron flux below '--^10 MeV is mainly due to sponta- 
neous flssion (the most important source being ^ssy ^ (a,n) processes due 
to interactions of a's from natural emitters with light target nuclei in the rock. 
Neutrons with energies above ~1G MeV are produced by nuclear reactions in- 
duced by cosmic ray muons in the rock and their flux can extend up to a few 
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Figure 2; Gamma ray spectrum measured at LNGS (Hall A) with a small Ge diode. 



Table 1: Gamma-ray flux (7/scm'^ ) in the underground Hall A of LNGS. 



Energy interval 


Gamma flux 


[keV] 


[7/s cm^ ] 


- 500 


0.51 


500 - 1000 


0.12 


1000 - 2000 


0.08 


2000 - 3000 


0.015 
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Table 2: Neutron flux measurements at the Gran Sasso laboratory [13. [l^. 



Energy interval 


Neutron fiux[10-'5 n/scni^] 


[MeV] 


Ref. [13j 


Ref. [14] 


10-3 - 2.5 




0.54 ± 0.01 


1 - 2.5 


0.14 ± 0.12 




2.5 - 5 


0.13 ± 0.04 


0.27 ± 0.14 


5 - 10 


0.15 ± 0.04 


0.05 ± 0.01 


10 - 15 


(0.4 ± 0.4) 10-3 


(0.6 ± 0.2) 10-3 


15 - 25 




(0.5 ± 0.3) 10-6 



GeV. Results obtained by neutron flux measurements at LNGS fl3, IS] are re- 
ported in Tab. [2] Simulations for both radioactivity induced and muon induced 
neutrons are found in literature [l5 , 17 1 . In the case of radioactivity induced 



neutrons the comparison with experimental results is rather good, while in the 
case of muon induced neutrons this comparison is not possible since this compo- 
nent dominates the environmental neutron flux only above ^10 MeV, a region 
where poor experimental data is available. This flux can be however computed 
on the basis of Monte Carlo simulations. Exis ting simulations performed to 
evaluate the muon induced neutrons at LNGS (isl [l7| yield an integral flux 
above 10 MeV that is three orders of magnitude lower than the radioactivity in- 
duced one, that amounts to ^A-10~^ n/scm^ below 10 MeV. Also muon induced 
cascades or showers in experimental materials (especially high A materials for 
gamma-ray shielding such as lead) can be an important source of hard neutrons. 

Neutrons are absorbed after moderation, this step being generally obtained 
by using several cms layer of hydrogen rich compounds (e.g. water and polyethy- 
lene). Neutron absorption is instead based on high neutron capture cross section 
isotopes (B and Li are among the most used). 



3.3. Muons 

The cosmic muon flux at LNGS is reduced by about six orders of magnitude 
with respect to the outside flux by the 3600 m.w.e. shield, provided by the 
mountain itself The average muon energy in the underground site is 270 
GeV [3 with an integrated flux of -(3.2 ± 0.2)-10-^ /i/scm^ 19] and average 
zenith angle {0) — 35°. The distribution of the azimuth angle ((> follows the 
proflle of the mountain. The differential flux and zenith angle dependence have 
been measured by LVD 3] and MACRO 20]. 

Generally, direct interactions of muons in the detector are not particularly dan- 
gerous. A minimum ionizing particle traversing a CUORE crystal would cause 
an energy deposition of a few tens of MeV, which is well above the Q- value of 
f5P[Qv) . Moreover the contribution to the background of direct muon interac- 
tions as well as of the muon generated electromagnetic showers, can be reduced 
by means of an anti-coincidence cut between different detectors. On the other 
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hand, gamma and neutron background arising from muon interactions in the 
setup are very dangerous. The presence of high density and high A shields 



makes the situation even worse (the neutron yield is proportional to A'^-^ 21 



even if simulations of muon produced neutrons in such materials appear to be 
subject to a large uncertainty. 

Muons, being deep penetrating particles, cannot be removed by a shield. 
The only way to get rid of this background, in those experiments in which it 
could be dangerous, is by the use of a muon veto to reject any muon induced 
interaction. 



4. Geant4 simulation 

Monte Carlo simulations play a primary role in background identification 
and therefore, in its reduction. Understanding the background event rate in 
an experiment relies on accurate Monte Carlo modelling of the experimental 
apparatus and Physics processes leading to all types of background from all 
possible sources. 

In this work, we present the results obtained with the version 4.8.2 of the 



GEANT4 [22[ simulation tool. GEANT4 is widely used among the particle 
physics community and its reliability has been established in a number of vali- 
dation studies. 

The LBE (Low Background Experiments) physics list has been used. It is 
suitable to handle both low energy interactions and hadronic interactions. It 
contains the Standard and Low Energy electromagnetic processes as well as a 
complete description of hadronic interactions at low and high energy. The High 
Precision (HP) neutron models are used to model capture, elastic scattering and 
inelastic scattering of neutrons from thermal energies up to 20 MeV. Gamma, 
electron and positron production thresholds for secondaries have been set to 5 
cm for lead and to 500/.tm for copper and iron. 

The CUORE geometry used in the simulation is shown in Fig. |31 The setup 
of the external shield has been simplified with respect to the real CUORE de- 
sign, described in section [2. II It consists of an octagonal structure surrounding 
the cryostat, with a 20 cm neutron shield as the more external layer and a 25 cm 
octagonal Pb shield inside. The 988 detectors are arranged in 19 vertical towers, 
each tower consisting of 13 layers of four cubes each, as in the real CUORE de- 
sign. The copper frame structure used to assemble the detector towers and the 
wire and source trays are also implemented. The small detector details (pins, 
thermistors, gold wires, teflon spacers, glue spots, etc.) are not implemented, 
due to their negligible influence on the results concerning environmental contri- 
bution to the experimental background (a few tens of grams total mass against 
the hundreds tons mass of the shielding). 

The output of the Monte Carlo simulation is converted in detector output 
by means of an additional code, G2TAS [l3|. In order to take into account 
the features of the detector response, namely the slow pulse development and 
the dead-time, two events depositing energy in the same detector are considered 
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Figure 3: 3D view of the geometry implemented in the simulation: external octagonal PET 
shield (left) and internal shields and detector (right). 



as a single event if the time interval between them is less than an integration 
time of 10 ms. A decay rate can be specified for each source, than all pile-up 
events (two events in the same bolometer within the acquisition window of 4 
s) are rejected. Coincidence patterns between the detectors (based on event 
multiplicity, evaluated both in the entire detector array and for the nearest 
neighbors) are invoked when computing spectra and scatter plots, as it is for 
experimental data. An energy threshold of 50 keV is assumed. The detector 
resolution is taken into account introducing a gaussian smearing with 5 keV 
FWHM. The ROI has been chosen as a symmetric region 40 keV wide around 
the nominal Q-value of the l3P{0v) of ^^OTe ((2527.517 ± 0.013 keV). In this 
region no gamma photoelectric peaks are expected (it has been chosen above 
the ^°Co sum gamma line at 2505 keV and below the 2614.7 ^°^T\ line). 



5. Simulation results 

A full Monte Carlo simulation for muon induced background should be per- 
formed generating muons in the rock and propagating them and all secondary 
particles produced in the rock and in the setup materials to the detector. In 
the present work, we adopted a different approach: 

• Muons have been generated and propagated from the cavern boundary; 

• Neutrons induced by muon interaction in the rock have been separately 
generated and propagated also from the cavern boundary; 

• High energy gamma rays induced by muons in the rock have not been 
simulated since the evaluated flux is much smaller than the one due to 



rock radioactivity 11 1 
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The chosen approach leads to an overestimate of the background due to the loss 
of correlation between neutrons (and other secondaries) and the original muons. 
This can cause a multiple site event to be identified as a single site one, thus 
reducing the anti-coincidence cut efficiency. 

The contribution of the environmental gammas to the background in the ROI 
has been evaluated by simulating primary gamma events at 2614.7 keV (indeed 
the dominant gamma contribution to the background comes from the ^osr^j 
gamma line, the strongest natural gamma line above the Q- value). According 
to the spectrum measured in the Hall A of LNGS (see section I3.ip the flux at 
this energy is 9-10~^ 7/scm^ . The primary gamma events have been uniformly 
generated on a spherical surface of 500 cm radius, surrounding the most external 
shield. The background value in the ROI is obtained by multiplying the integral 
of the Monte Carlo spectrum in that region by the normalization factor between 
the number of generated events and the actual flux of 2615 keV gammas at 
LNGS. 

The results are reported in the first row of Tab. [3] where the total, the global 
anti-coincidence and the near channel anti-coincidence contributions are shown. 
No anti-coincidence cut is applied in the first case, while an anti-coincidence cut 
between all the 988 detectors and between each detector and its 18 closest ones 
are applied in the second and in the third case, respectively. The probability for 
a gamma ray generated from outside the external shield to reach the detector is 
very low and, due to multiple Compton interactions, the energy released is low 
and below the I3(3{0i/) Q- value. No event has been recorded in the Monte Carlo 
in the ROI, so only an upper limit has been obtained by considering Poisson 
statistics, well below the CUORE background goal. 

The neutron spectrum has been simulated according to Fig. 21 The energy 
interval between 25 meV (thermal neutrons) and 1.5 GeV has been divided in 22 
logarithmic energy bins. Flux values are taken according to existing experimen- 
tal dat a [isl. 14 [ below 10 MeV and to simulations of muon induced neutrons at 
LNGS [ii Q above 10 MeV. In order to be conservative, for each energy bin 
the highest flux among different measurements or different simulations has been 
used. The total neutron flux, thus extracted, is found to be: (3.7 ± 0.2) • 10~^ 
n/s cm^ . The results are summarized in the second row of Tab. [3] Even the total 
background is well below the CUORE requirement. Furthermore, a reduction 
factor greater than 1 order of magnitude could be obtained by applying a global 
anti-coincidence cut and a reduction factor of ~5 by applying a near channel 
anti-coincidence cut. In Fig. [6] a comparison between the neutron background 
spectrum before and after the global and near channel anti-coincidence cuts is 
shown. Simulations performed splitting the neutron spectrum in the energy bins 
of Tab. m show that the largest contribution is due to neutrons with energies up 
to 10 MeV, i.e. radioactivity induced neutrons. 

Muon induced background was modelled by generating muons on a spherical 
surface of 500 cm radius (similar to gamma rays and neutrons) and transporting 
these muons along with all secondary particles produced in muon interaction in 
the detector materials, through the whole setup. All energy depositions in the 
crystals were recorded. The muon spectrum angular distribution is parameter- 
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Figure 4: Simulated neutron flux density as a function of the neutron kinetic energy. Values are 
taken accordingly to experimental data ll^ [l^ (below 10 MeV) and to existent simulations 
of muon induced neutrons at LNGS (above 10 MeV). 



ized on the MACRO data (see sectio nl3.3p . The energy spectrum is generated at 
the ground level (GL) accordin g to [23j, that in the high energy limit {Eol > 1 
TeV) can be approximated to [18|: 



dN 



0-14 ■ E^l' 



0.054 



dEoL ■ s crn^ sr GeV v 1 + a Eql cos 9 1 + /? Eql cos 6 



(3) 



where a=l. 1/115 GeV and /3=1. 1/850 GeV. It is then extrapolated underground 
using the equation [23| : 



Eu = {Egl +e)-E- 



-bX 



(4) 



where b=0.4-10^^ g^'^ cm~^, e=540.0 GeV and X is the thickness of the rock 
traversed by the muon. X is computed for each given direction taking into 
account amount and shape of the overburden mountain according to MAGRO 
data. The simulated muon energy spectrum is shown in Fig. [5 The energy 
mean value is consistent with the measured one of 270 GeV jlSl, |24|. 

The muon energy has been chosen between 1 GeV and 10 TeV underground, 
including 99.2% of the total muon flux. The resulting total background from 
cosmic ray muons (third row of Tab. [3]) is higher than the GUORE goal. A 
reduction factor of almost two orders of magnitude is obtained by applying a 
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Figure 5: Simulated muon energy spectrum. 



global anti-coincidence cut and a reduction factor of ^10 by applying the near 
channel anti-coincidence (Fig. [6]). The CUORE background requirement could 
thus be reached without the necessity of a muon veto system. 

The errors quoted in Tab. [3] and for the background results are statistical 
only. The neutron flux above 10 MeV, due to cosmic muon interactions in the 
rock surrounding the laboratory, has a systematic error of about a factor 2. 
This value takes into account both the discrepancies between the Monte Carlo 
calculations performed with different codes (i.e. Geant4 and Fluka) [2E\ and 
between the calculations and the available neutron yield data for light targets 
and for heavy targets like lead [21]. Below 10 MeV the systematic uncertainty on 
the measured neutron flux is evaluated to be 6.5% for [lj| and 31.6% for [l^] and 
the two measures have discrepancies up to a factor 3. Also with the conservative 
approach adopted (see above) the event rate is much below the CUORE target. 
An additional source of systematics in the background evaluation is due to 
neutron propagation and detection in the simulation code, mainly related to the 
precision of the cross section tables used to model capture, elastic and inelastic 
scattering. According to [llj it has been assumed to be 20%. 

The systematic error in the muon flux is 8% [l9| . Uncertainties in the muon 
induced neutron yield, neutron transport and detection must be also taken into 
account as a source of systematic error together with the uncertainties in the 
models used in Geant4 to evaluate the electromagnetic component. The former 
has a total error of 45%, the latter of 5% [11|. 
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Table 3: Background in the ROI obtained by simulations of the evaluated gamma, neutron and 
muon fluxes at LNGS. Quoted errors are statistical only. Units are: counts ■ 10^'^ /kg keV yr. 
Limits are given at 90% C.L. 



Source 


Statistics 


Total 


Anti-coinc. 


Anti-coinc. 




[years] 




(Global) 


(Near chan.) 


Gamma 


0.1 


0.390 


0.390 


0.390 


Neutron 


7.88 


0.270 ±0.022 


(8.56 ± 6.06) X 10-3 


0.0642 ± 0.0442 


Muon 


4.05 


17.3 ±0.3 


0.104 ±0.022 


1.850 ±0.049 




Figure 6: Comparison of the Monte Carlo spectra obtained by simulating neutron flux (left) 
and muon flux (right) at LNGS. In each plot the spectrum obtained without any anti- 
coincidence cut (green) is shown together with the ones obtained by applying the global 
anti-coincidence (blue) and the near channel anti-coincidence (red) cut respectively. On the 
right side of each plot the axis reporting the expected counts/kg keV yr in the CUORE detector 
is shown. 



The uncertainties in the environmental gamma flux at LNGS originate from 
the measurement statistics and from the detector efficiency evaluation, per- 
formed with Monte Carlo simulations. The gamma peaks are measured with 
uncertainties of less then 8%, while the simulation uncertainties in the Geant4 
electromagnetic tracking have been evaluated to be at a few percent level jll|. 
Gamma transport, interaction and detection uncertainties in the Monte Carlo 
code used to extrapolate the CUORE background are an additional source of 
systematics at a few percent level. 



6. Conclusions 

The contribution of environmental and cosmic radioactivity at LNGS (i.e. 
neutrons, gammas, and cosmic muons) to the expected background in the 
CUORE ROI has been evaluated by means of simulations performed with a 
GEANT4 based code. Even without the use of a muon veto, the chosen shield 
configuration guarantees that the overall contribution of these sources is less 
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than < 10 ^counts /kg keV yr by operating detectors in near channel anti-coincidence, 
a value that satisfies the CUORE background goal. 
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